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Abstract 

We investigate Gibbs measures relative to Brownian motion in the case when 
the interaction energy is given by a double stochastic integral. In the case when 
the double stochastic integral is originating from the Pauli-Fierz model in nonrel- 
ativistic quantum electrodynamics, we prove the existence of its infinite volume 
limit. 



1 Preliminaries 

1.1 Gibbs measures relative to Brownian motions 

Gibbs measures relative to Brownian motion appeared in [22], where they have been 
introduced to study a particle system linearly coupled to a scalar quantum field. A 
systematic study of such measures has been started from [B], where by making use of 
this measure the spectrum of the so-called Nelson model is investigated. Since then 
there has been growing activity and interest in the study of various types of these 
measures H [121 [13 E] • 

One way to understand Gibbs measures relative to Brownian motion is to view 
them as the limit of a one-dimensional chain of unbounded interacting spins, with the 
distance between the spins going to zero. As a simple example, which will be instructive 
in what follows, let us take R'^ for the spin space, and fix a (finite or infinite) a priori 
measure uq on R'^ as well as smooth, bounded functions V : R*^ — * R and W : R x R'^ — > R. 
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Double stochastic integrals 



On the lattice eZ fl [— T, T] with spacing e and n = 2T/e sites, we define the measure 
z/^ through 

i/*(dx_„ . . . dxn) = ^ n i^o(dx,)e-^^^^("')-^^»("'+^-">)'+^'^-^-*(^-*'"^-"'). (1.1) 



|j|<n 



Here Z,^w normahzes u'^ to a probabihty measure, and for finite e, u'^ is just a chain 
of interacting spins. However, the scahng becomes very important when e ^ 0. Then, 
formally each spin configuration {xi)\i\<n € eZ fl [— T, T] becomes a function x(-) on 
[— T, T] , and the particular scaling of the quadratic term above gives rise to the term 

limg^o ^ Xli i^i+i ~ ^li^ 1^^ = J^j,{dx{s) /dsyds. It is this term that prevents the mea- 
sure z/* from being concentrated on more and more rough functions when e —>■ 0, 
ensuring continuity of x{t) in the limit. Indeed, when uq is chosen as the Lebesgue 
measure on R'^, it is not difficult to show that part of the normalization along with the 
quadratic term give converge to Wiener measure W, so that in the limit, e 0, we 
obtain 



Zy« 

Here {Bt)t>Q is now a Brownian motion, hence we call (11.20 Gibbs measures relative 
to Brownian motion. Indeed, the measure appearing in [22] is of the above type, and 
most of the subsequent works cited above have been concerned with measures of the 
form ffT:2|) . 

In this paper we study another type of Gibbs measures, arising from a very similar 
discrete spin system. Namely, let us now define 

i.^-^(da:_„...dx„) 

Yl Mdxi) exp I -e ^ V(xi) - - 



\i\<n 



+ ^(xi+i - Xi) ■ Wa/(j - ^, Xj - Xi){xj+i - Xj) . (1.3) 



Now Wm is a dxd matrix, but otherwise the expression looks very similar to (11. ip . The 
crucial point is, however, that now the scaling of the term involving Wm is different. 
The e"^ which ensured convergence to a double Riemann integral is gone by sandwiching 
e^Wm between (xj+i — Xi)/e and (x^+i — Xj)/e, and replaced by the increments of the 
spins themselves. Since these increments will eventually converge to Brownian motion 
increments, as discussed above, they are of order i/e, so the scaling is indeed different. 
So after taking the limit e — > 0, we informally obtain 

V^'''{dB) = I-T'^(Bs)ds-J^^ f^^dB,-WM{Bt-Bs,t-s)dBt^yY. (1.4) 

As a consequence, taking the limit e — > yields a double stochastic integral in place of 
the double Riemann integral (II. 2p . 
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1.2 Definition of double stochastic integrals 

From now on we assume that d = 3 and specify the pair potential Wm = W = 
W{X,t) = (l^;..(X,t))i<^,,<3 given by 

^^'^^^'^ ■= I (1-5) 
where S-^{k) = i5^uik))i<^,u<3 is given by 

^tM ■■= V - (1-6) 



Measures of the type (11.41) with pair potential (II. 5p appear in the study of the so-called 
non-relativistic quantum electrodynamics, and have been introduced on a formal level 
in [T0l[T^[2^ . However we notice that there are some difficulties in the expression (11.41) : 
For t > s, the integrand is not adapted to the natural filtration = a{Br]r < T), 
so as a stochastic integral or any of its obvious transformations the double stochastic 
integral such as (ll.4p does not make sense. So the right-hand side of (11.41) is just an 
informal symbol. 

In [T6l Definition 4.1] and [HI (4.2)], however, the firm mathematical definition of 
(ll.4p has been given through a Gaussian random process associated with an Euclidean 
quantum field. We outline it below. A Gaussian random process £^^{f) labeled by 
/ G ©^L^(R^+^) on some probability space (<^e, Sg, ixe) is introduced, which has mean 
zero and covariance ¥.^j^[£^{f)£^{g)] = q{f,g) given by 

QiL 9)--=\j TiKko) ■ S^mik ko)ci'kdko (1.7) 
for f,g& ©^L^(R^+^), where ' denotes the Fourier transformation. Let 

Kt = (BU f JsVi- - B,)dB^ (1.8) 
Jo 

be the ©'^L^(R^+^)-valued stochastic integral defined in the similar way as standard 
stochastic integrals, where js : L^(R^) -h> L^(M'^+^) denotes the isometry satisfying 

ijsfJt9)L^R3+i) = (/,e~l*"^l'^^)i2(R3). (1.9) 
See Subsection 3.1 for the details. 

Definition 1.1 Let W be the pair potential defined in ( fi.5]) . The dounle stochastic 
integral is defined by 

[ [ dBs-W{Bt- B,,t-s)dBt:=qiKt,Kt). (1.10) 
Jo Jo 

We would like to express (II.IOD as an iterated stochastic integral in this paper. 
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1.3 Main results 

Let us define the Wiener measure W on X := C{R; R^), cf. also [231 P- 39, Remark 1.]. 
Let Hq = — (1/2)A. Suppose that /i,...,/n-i € L°°(R^) with compact support. Then 
there exists a measure W on X such that 

= / MK)MBt,)---fn{BtJdW. (1.11) 
Jx 

A path with respect to this measure is denoted by Bt{w) = w{t) for w E X. Note 
that Wiener measure is not a probabihty measure, indeed it has infinite mass. If 
denotes the measure of standard Brownian motion starting from x G at time to, 
then 

/ fo{Bj---UB,JdW= [ dx [ MB,,)---UB,JdP^'". 

Jx Jr3 JC([to,oo);K3) 

Let ip G L'^(R^) be a nonnegative function and we fix it throughout this paper. In the 
case of f ll.2l) . the existence of the weak limit of the measure on X, 

as T — > oo has been investigated for various kinds of V and W, and the limiting measure, 
z/^, proved to be useful to study the ground state ipg of some particle system linearly 
coupled to a scalar quantum field. Namely for a suitable operator O, we can express the 
expectation {ipg, Oipg) as fodi^^ with some integrand fo- So, beyond the existence 
of a measure of the form (11.41) . one is interested in the limit as T — > oo, at least along 
a subsequence. In other words, one would like to prove the tightness of the family of 
measures (11.41) . This is by no means an easy task, given that there are very few good 
general estimates on single stochastic integrals, let alone double integrals. 

The purpose of our present paper is to point out that there is at least one special 
case where there is a comparatively easy way to construct both the finite volume Gibbs 
measure and the infinite volume limit, namely the case when Wm = W = (W^iy)i<^^i,<3 
is given in (11.51) . Fortunately, this special case is the one that motivated the whole 
theory of Gibbs measures with double stochastic integrals. The main results in this 
paper are 

(1) we give an iterated stochastic integral expression of (11.101) : 

(2) we show the tightness of the family of measures 

BtM B_T)e- nBs)ds-a^ dB,-W{B,-B,,t-s)dB,^y^ 

for a general class of V including the Coulomb potential V{x) = —l/\x\, and arbitrary 
values of coupling constant a G K. 
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There has been recent progress both of the above topics: M. Gubinelh and J. 
Lorinczi [13] employ the concepts of stochastic currents rough paths in order to define 
f ll.4p rigorously for finite volume, and use a cluster expansion in order to construct the 
infinite volume limit. While these are impressive results, the techniques used are rather 
advanced, and the use of cluster expansion comes with strong assumptions on single 
site potentials V and coupling constants. The advantage of our methods is that we 
can avoid some restrictions needed in [13]; in particular we need not restrict to single 
site potentials that grow faster than quadratically at infinity, and we need no small 
coupling constant in front of the double stochastic integral. In particular our results 
include the Coulomb potential which is the most reasonable single site potential. On 
the other hand, of course the range of potentials W that is treated in [13] is much 
greater than ours. 

The paper is organized as follows: In Section 2 we will construct the finite vol- 
ume Gibbs measure as the marginal of a measure with single stochastic integral on 
a larger state space. This construction is well known [25], but has not been carried 
out rigorously so far. In Section 3, we rely on the detailed results available about the 
Pauli-Fierz model [Sni [H] in order to show that our family of Gibbs measures is tight, 
giving the existence of an infinite volume measure. While we expect that the general 
method of enlarging the state space should allow us to define and prove infinite volume 
limits for many more models than just Pauli-Fierz, this is not all straightforward. We 
will comment on this issue at the end of Section [3l 

2 Iterated expression of finite volume measures 

In this section we will specify the measure fiT that we are working with, and identify 
it as the marginal of another measure on a larger state space. Let us start by 
introducing an infinite dimensional Ornstein-Uhlenbeck process which will serve as the 
reference measure for the auxiliary degrees of freedom. Put 

uj{k) = ^\k\^ + m^ (2.1) 

for m > 0, and let Xs{f) be the Gaussian random process on a probability space 
(^, E, Q) labeled by measurable function / = (/i, /2, /s) with mean zero and covariance 
given by 

^G[Xsif)X,ig)] = J d'k^^e-^^'^^'-^W) ■ 5Hk)g{k). (2.2) 

Here / denotes the Fourier transform of / and we assume that /^/ ^/uJ, g^/ ^JZ) E L^(M^), 
/i, z/ = 1, 2, 3. 

Remark 2.1 Let Ys{f) be the Gaussian random process on {^e, ^e^IJ'e) defined by 

Ysif) ■■= ^^(j.(/7v^)"). (2.3) 
Then Ys{f ) is mean zero and its covariance is 

E,,[Umif)] = Eg[X,if)X,ig)]. (2.4) 
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Hence Ys{f) and Xs{f) are isomorphic as Gaussian random processes. 

We will now couple Q to the Wiener measure W. For this we use a coupling function 
ip with the assumption below: 

Assumption (A): 

(1) 0{k) = 0{-k) = 0{k) and ^JZj(p,(pl i^j G 

(2) (p is rotation invariant, i.e., 0{Rk) = 0{k) for all R G 0(3). 
Let us now define the quantity 

J[o,T](X) := / X,((p(- - B,)) ■ dB,. (2.5) 
The proper definition of J[o,t] reads 

n 

JlO,T]{X) := lim y X(^j_i)T/n{fi- - B(j^i)T/n){BjT/n " -B(j_l)T/n.) ) , (2.6) 

where the right hand side strongly converges in L2(X X ^; ^ ® P^""). This is proved 
by showing that the right-hand side of (12. 6p is Cauchy by making use of (12. 2p . In the 
same way, we can define 

J[_T,T](X) := - Bs)) ■ dB,, (2.7) 

where (i?t)teE denotes the 3-dimensional Brownian motion on the whole time line K. 
The coupling between the Gaussian process and Brownian motion is given by the 
measure u on X x ^ with 

dz/T = ^ exp (^ia J X,(^(- - Bs)) ■ dB^^ ij{B^T)i^{BT)dW ® dQ, (2.8) 

where ip G L^(M^) is an arbitrary nonnegative function, Zt the normalizing constant, 
and a is a coupling constant. In order to guarantee that the density in (12.81) is integrable 
with respect to W, we chose the boundary function t/.' to be of rapid decrease at infinity. 

We are now in the position to define our finite volume Gibbs measure. We will 
introduce an on-site potential V which we take Kato-decomposable [7], i.e. we require 
that the negative part V- is in the Kato class while the positive part V"+ is the locally 
Kato class [23]. This ensures e.g. that 

sup Ep=:,0 



exp 



V{Bs) ds 



< oo. 



(2.9) 



Double stochastic integrals 



7 



Definition 2.2 Let V : M. be Kato- decomposable and a G M a coupling constant. 

Then the measure fi^ on X is defined through 



(2.10) 



We want to show that the measure /x^ we just defined is a Gibbs measure with 
double stochastic integral as given in Section [H The key to doing this is the fact that 
we will be actually able to calculate the Gaussian integral exp(iJ[_2-2-](X))d^(X), 
and thus are left with an expression involving Brownian motion paths only. In doing 
so, we will set a = 1 for a simpler notation. 

Let us give the heuristic presentation first. By the standard formula we have 

E^[eiJ[-T,T,] = exp ('-^Eg[jf_^,^]]') (2.11) 



and formally, by Remark 12. H we have 

Eg[J[lT.T]] '°=" If / ■ W{Bt -B,,t- s)dBt, (2.12) 

2 J_J' J_rp 

where W is given in (11. 5p . As it stands, there are problems with the right-hand side of 
formal expression (I2.12p . mainly because the integrand is not adapted. The resolution 
is to use symmetry of W and break up the integral into two parts, one where s < t and 
one where s > t, which are then proper iterated Ito integrals. This leaves the diagonal 
part, which gives a non- vanishing contribution by the unbounded variation of Bt. 
We define the iterated stochastic integral St by 

Sj. := j d^k^-^^ e* »-dB. • j' e-"<'>l'-"e-"--a*-^(t)dB, + 

+1 " (2.13) 

3 J 2uo{k) ^ ' 

St is the well-defined expression that will replace (12.101) . The above line of reasoning 
and (I2.13P are not new [25], except that (I2.13p is usually not written out but instead just 
referred to as the double stochastic integral with the diagonal removed. Nevertheless, 
(I2.13P can be considered as known. However, the derivation above is mathematically 
not rigorous, since the ill-defined expression (12.120 appears along the way. To avoid 
this, one has to derive (I2.13p directly from Eg[e'''[-^'^i]. This is what we do in the next 
theorem. 

Theorem 2.3 For almost every w E X , we have 

Eg[e'-^[-^'^i] = e"^^. (2.14) 
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Proof: Let us replace the time interval [— T, T] with [0, T] for notational convenience. 
We employ (12. 6p and use dominated convergence to get 



lim exp 



J=l 



where we set SBj = Bjx/n — -B(j_i)T/n and = (j — l)T/n, j = 1, A^. Now 



Er 



J2xaM--Ba,))-sb, 



.i=i 



2V / d^fc^^l^e-^^-^^^e^'^^-^ V (2.15) 



1=1 



^3^M^5^(A;))55^, (2.16) 



For the diagonal term in the last line above we note that 



2 f\mi 

2uj{k) 



d^k 



by the rotation invariance oi0. Now as n ^ oo, Yl^=i I'^-^iP ~^ ^5 while Yl]=i l^^jl"^ 
T/3 for all /i = 1, 2, 3, for almost every w ^ X. Thus for almost every G X, we find 



i=i 



2uj{k) 



2T 



2uj{k) ' 



For the off-diagonal term, we start by noting that by the definition of the Ito integral 
for locally bounded functions f , g : R x ^ R, we can see that 



]EpO,0 







2" 


/ ds 


f{s,B,) rg{r,Br)dB^ 




Jo 


Jo 





< 00. 



Hence the stochastic integral of p^{s) = /(s, Bg) g{r, Br)dBI^ exists for all yU = 1,2, 3, 
and it holds that for each k ER^, 

lim V /(A„ Ba^)8Bj ■ 5^{k) / g{r, Br)dB, 
"^^^ ~t\ Jo 

= £ f{s, B,)dB, ■ 6^{k) (^^' g{r, B,) dB, 



(2.17) 
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strongly in L'^{X; P^^). By the independence of Brownian increments and the fact 
that EpO.o = 1/n, KpO,o[SBj] = 0, we can estimate the L'^{X; )-difference of 
fl2.17p and the off-diagonal term: 

2 



EpO.O 



2 1 

3 n 



3=1 V'> 1=1 



2 -, 3 n 



/ g{r, Br) dB'; - ^ g{Ai, Ba,)5B^i 



(2.18) 



Then the right-hand side above converges to zero as n — > oo, and together with fl2.17p 
we find 

n j 

hm ^ /(A„ Ba,)SB, ■ 6^{k) J2 9{Ai, Ba,)SBi 



i=i 1=1 

f{s,B,)dBs-6^{k) [ g{r,Br)dB, 
Jo 



(2.19) 



in L2(X;P^°). By putting f{t,x) = e'k-x^-u,{k)t g(^^^^^ ^ ^-ik-x^u,{k)t -^^ ^rm . we 
can see that fl2.15p converges to the off-diagonal part of 5^. Then the proof is finished, 
qed 

Remark 2.4 It is interesting that we know that [e""^^! = |Eg[e''^'-'^''^i] | < 1 almost 
surely. This is not obvious from the iterated integral representation e""^"^. 

Let us summarize: 

Proposition 2.5 Let fi}f be the measure on X from Definition \2.2[ Then 
d/i^ = -l^/,(p_^)^/,(Py)e-"'^^e-^-T^(^«)<^^ dW, 

where St is defined by St with the diagonal part removed: 



St := / d^A; 



f d^A;^^ f e^'^-^^dP, ■ r e--('=)(^-'-)e-i^-^'-5^(A;)dP,. 
J 2uj{k) J_T J-t 



Or 



where 



St '■- 



Z{s, w) ■ dBg 



Z{s, w) 



dBr 



2uj{k) 



5^{k)e-^'-''^''^^'^e-'^-^^^-^^^d?k I . 
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Remark 2.6 In Proposition 12. 5[ the diagonal term |^ / d^k ^^^^H is absorbed in the 



normahzation constant, since it does not depend on the Brownian path B. Moreover 
from Remark 12.41 it follows that 



exp(-..)|.e.p(I/.3.M|) 



Thus St may be written symbolically as 

St = \ j dBs- W{Bt -Bs,t- s)dBt. 

3 The infinite volume limit 

3.1 Tightness and the Pauh-Fierz model 

The idea of the proof of the infinite volume limit we are about to give is not straight- 
forward. We will show that it follows from showing that the bottom of the spectrum 
of a self-adjoint operator is eigenvalue. Actually, in the case of pair potential W under 
consideration, associated self-adjoint operator is realized as the Pauli-Fierz Hamilto- 
nian B. in the non-relativistic quantum electrodynamics. Fortunately it is established 
that B has the unique ground state for not only confining external potential V ^ e.g., 
= but also the Coulomb Vi^x) = —l/\x\, which is the most important case. 

Let us begin with defining the Pauli-Fierz Hamiltonian with form factor a.s a. self- 
adjoint operator on some Hilbert space and we will review the functional integral 
representation of the Co semigroup e~^^ . 

Let ^ := 0r=o[(S)"^^(]^^ X {1>2})] be the Boson Fock space. The state space of 
one electron minimally coupled with the photon (bose) field is given by 

:= L'^{R^)(g)^. (3.1) 

We denote the formal kernels of the annihilation operator and the creation operator 
on ^ by a{k,j) and a*{k,j), respectively, which satisfy the canonical commutation 
relations: 

[aik,j),a*{k',j')] = 5ik-k%,,, [aik,j),a{k',j')]=0 = [a*{k,j),a*{k',j')]. (3.2) 
The free Hamiltonian in ^ is defined by 

Br.= J2 I ^{k)a*{k,])a{k,])d?k. (3.3) 
i=i,2 

Here dispersion relation uj is given by 02. II) . Let us fix a function satisfying Assump- 
tion (A) The quantized radiation field A = [Ai, A2, A3) with form factor is defined 
by A^ := f®^ A^{x)d^x, where we used the isomorphism = J^^ L^(R^) da; and 



j=l,2 
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The vectors e{k,j), j = 1, 2, are the polarization vectors. They satisfy e{k, i) ■ e{k,j) = 
6ij and k ■ e{k,j) = 0. Note that 

J2e,{k,j)eAk,j) = Si^{k). (3.5) 

3=1,2 

(13. 5p is of course independent of the choice of polarization vectors and k ■ e{k,j) = 
yields that 

3 

5^V.,A^(x) = 0. (3.6) 

The Pauli-Fierz Hamiltonian H{0) is defined by 

H{0) ■.= ^{-tV ®l-aAf + l®Hi, (3.7) 

where a G M denotes coupling constant. It is established in [19] that H{0) is self- 
adjoint on D{—A) n D[Hf) and bounded from below. Moreover H{0) is essentially 
self-adjoint on any core of— (1/2)A(8)1 + 1(8) iff. We now introduce a class of external 
potentials : ^ M that we can add to Hq. 

Definition 3.1 V e K if and only ifV = V+-V^ such that V± > 0, V+ e Ll^iR"^) 
and relatively form bounded with respect to —(1/2) A with bound strictly smaller 
than one. 

Let V & K. Then we define H as 

H := H{0) + V+0 1 - V-(g)l, (3.8) 

where ± denotes the quadratic form sum. To see the weak convergence of yU.^, we 
introduce the assumption below. 

Assumption (GS): There exists a ground state (fg of H. 
Example 3.2 Let 

Vix) = ~^ + Uix), (3.9) 
\x\ 

where C > is a constant, and U = U+ — G Lioc(lR'^) such that U± > 0, 
infa-gRS U{x) > — oo, t/_ is compactly supported, and —(1/2) A + U has a ground 
state (p > with ground state energy — cq < such that |0(x)| < ^e"'^''/'^ with some 
constant 7 > 0. Then the ground state of H exists for arbitrary values of a. See [H 
p. 8] and [HITT]. Typical examples are 

C 

^Coulomb ( 2^) = "TT' 

\x\ 

Confining (•^) l-^l ) ^ 1)2,... 
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To construct the functional integral representation of e *^ we introduce some proba- 
bilistic notation which was already mentioned in Section 1. Let {■^^{f)}j^^3i^2ij^3+iy 
denote the Gaussian random process labeled by / G ©^L^(]R^"''^) on some probabil- 
ity space {^E,'^E, ^J'e) with mean zero and covariance given by K^j^[s^^{f)£/^{g)] = 
q{f,g), where g(-, ■) is defined in (11 ■71) . We define the isometry jg : L^(M^) L^(m^+-'^) 
^ ''"'y^/^)^yuJ{k)/{uJ{ky + |A;oP)/(A;) which satisfies ([L9]). The cru- 



by jsfik, ko) := (e 
cial identity linking the Pauli-Fierz model to Gibbs measures is 

Proposition 3.3 

(1) For arbitrary f G L'^{R^) with f > but f ^ 0, it follows that 

((^g, / ® > 

(2) Let /i, G and i' G L'^{R?). For -T = to < ti < 
Euclidean n-point green function is expressed as 



(3.10) 

< tn = T, the 



{ip ® n, e"(*i^*»)^(/i ® 1) ■ ■ ■ {fn-1 ® l)e-(*"-*"-i)^^ ® n) 



E 



Proof: See [H] for (1). In [151 [20] it is estabhshed that 

® n, e-(*^-*«)^(/i ® 1) ■ ■ ■ (/„_i ® l)e-(*"-*"-^)^V^ ® f]),^ 



n-1 



U=i 



E- 



where 



MB 



1,2,3, 



-ia J^^ .s/f -dBs 



S G 



and A = Since = (?/) ® Q,e~^^^4j ® fi).^' and e 



Eg[e^"-^tO'^l 



E 



the lemma follows. 



qed 



Remark 3.4 Formally, (2) of Proposition 1X51 can be deduced from using the Feynman- 
Kac-Ito formula P [H [IS [201 |23l [25] but note that integrand depends on time s 
explicitly; although this formula would give the Stratonovitch integral jj o d-B^ = 
s^f ■ dBg — I /J V ■ jz/fds instead of the Ito integral s^^ ■ dBs above, the Coulomb 
gauge (13. 6p allows us to use the Ito integral instead, since V^; ■ j2^^(A(- — x)) = 0. 



By (I3.10p . we know that the the ground state, (fg, of H is unique if it exists and, in 
particular, (yPg, / ® fi)^ 7^ holds, then we can define the sequence converging to the 
normalized ground state ifg by 



-tH 



(/®l])||^ie-*^(/®fi). 
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Actually, by virtue of flS.lOp . we see that 

ipg = s - lim (fl. (3.11) 

One immediate and useful corollary of (13.111) and Proposition 13.31 is as follows. 
Corollary 3.5 Let p,pi,p2 e L°°(R'^). Then for t > s, 
lim E v[p{Bo)] = ((/5g, (p ® l)v?g),_^, 

1 — »-oo ^ 

lim E^v[p^iB,)p,iB,)] = (y,g, (pi ® l)e-(*-)^(p2 ® l)¥.g),^'e(*-)^(^), 
where E{H) = inf cr(i/) denotes the ground stare energy of H . 

In order to prove the main theorem, we show a more general formula than (2) of 
Proposition 13. 3[ Let 



{k)a*{k,j) + f^{-k)a{k,j))dk. 



Define the family of isometries J* : ^ — > L'^{^e), t E R, hj the second quantization 

of j„ namely J* :A(/i) ■ ■ ■ A{fn): =:^^(jt/) ■ ■ -^^(jjn): and J^fi = 1, where i^: 
denotes the Wick product of ^. 

Proposition 3.6 Let F,G E and /i, fn-i E L°°(M^). For S = to < ti < ■ ■ ■ < 
t =T 

{F, e"(*^-*°)^(/i ® 1) ■ ■ ■ (/„„i ® l)e"(*"-*"-)^G),^ 

HfjiB,^) e~fsyiBs)ds^^^ [3sF{Bs)e-'^^s<'-^B.3^G{BT) 



Proof: See qed 
We are now ready to state and prove the main theorem of this paper. 

Theorem 3.7 Suppose that Assumption (GS) and Ii2.9\) . Then there exists a subse- 
quence T' such that the weak limit of p^, as T' ^ oo exists. 

Proof: By the Prohorov theorem, it is enough to show two facts: 
(1) lim supp^d^op > A) = 0, 



A— ►oo 



(2) for arbitrary e > 0, lim sup p^ max \Bt — Bs\> e\ =0. 

<5J,0 X \ \t-B\<s 



-T<s,t<T 



Using Corollary 13.51 we have 

> A) = {ipl, (X{|xP>A} ® 
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where xd denotes the characteristic function on D. Using the fact that ip'^ ipg 
strongly as T — > cxo and ||X{|a;|2>A}V'g|l^ ^ as A ^ oo, we get (1). For (2), assume 
that |t — s| is sufficiently small. It is enough to show that 



E^v[\B,-Bsn<\t-srD 



(3.12) 



with some constant D independent of T. To apply Propositions 13.51 and 13. 6[ we have 
to truncate the process Bj^ as 

—a, B^{w) < —a, 



(Snaiw) : = 

and define the truncated multiplication operator /i^, = 1, 2, 3, by 
Note that 



a, B^{w) > a. 



-af{x), < -a, 
x^fix), \x^\ < a, 
af{x), Xf, > a. 



\ha{x) - ha{y)\ < \x - y\ x, ?/ G M^ 
for all a > 0. Since /i^ is bounded, we can see that 



(3.13) 



Bt)a — (-Bs)aP"] 



2n 

k 



3 2n 

EE 

u=l k=0 
3 2n p -, 

1^1^ k 

u=l k=0 - 
3 2n 



{K ® 1) 



2n-k +tH T 



.if 



EE 



u=l k=0 



2n 
k 



-l)'Ew 
X E 



KiBo)) (KiB, 



2n-k 



\ha{Bo)-K{B,^,)\'-e-iryiBs)^s 



X E 



'MB 



< E 



w 



I -Bo — Bt^ 



Joe-^V^^(5o)e-*"^--^="-^^=J*-.e+*^y.J(i?, 

1 1/2 



xE^ 



1/2 
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where we e^*^<y9g is well defined for t < T, and we used Corollary 13.51 in the second 
equality, Proposition 13.61 in the third equality, fl3.13p in the fifth inequality, Cy '■= 

sup^gjj3 IEp^_ [e~^-^ ^(^r-)drj ^ ^ g^j^j (j^^ ^Yie constant such that 

EpO,o[\B, - Btl'''] = C4n\t - s\^\ 

Since ||e~*"'^(/)g || e'^^^^^'WipgW and ||e"'"*"'^(y9g || e^^^^'^W^pgW as T oo, we have 
D := supj. y^Cy||v9^||^e(*-^)^(^) < oo. Then we have 

E^v[\{B,)^ - {B,)^n < D\t - sr 

uniformly in a. Since the left-hand side above monotonously increasing as a t oo, the 
monotone convergence theorem yields (13.121) . Thus (2) follows. qed 

Definition 3.8 Let V E K and suppose Assumption (GS). Then the weak limit of the 
measure u^, on X is denoted by u^. 

Using the functional integration of e~*^, we can show the Carmona type estimate 
[8], namely ipg is spatially localized as follows: if = then ||(^g(x)||^ < 

(j^Q-C2\x\" ^ g^j^^ jf V{x) = —l/\x\, then ||v5g(a;)||,^ < C^e''^'^^^^ for some constants Cj. 
We have a corollary. 

Corollary 3.9 Assume that \\Lpg(x)\\^ < Ce~^^^^~' for some positive constants C,c and 
7. Then 

[ e'='^°''z/^(dtf;) < 00. (3.14) 
Jx 

Proof: Let p^(x) = I ^ ' Then (<^g, (p^ ® l)<^g)jr = JxpUBo)fJ'L 

follows. By the limiting arguments as m — > 00, we have (13.141) . qed 



3.2 Concluding remarks 

In this paper we have given one example where we can both make sense of the double 
stochastic integral and obtain the infinite volume Gibbs measure by coupling Brownian 
motion to an auxiliary Gaussian measure. The drawback of this particular example is 
that the Gaussian space is infinite dimensional, and the associated Hamiltonian along 
with the existence of its ground state is non-trivial, and so we have to rely on a lot of 
technology. It is conceivable that the same method should work in a much easier case, 
namely when the auxiliary Gaussian process is just the stationary one-dimensional 
(or n-dimensional) Ornstein-Uhlenbeck process. However, when trying this approach 
one notices that on the way we used a lot of special features of the Pauli-Fierz model 
and its associated functional integral: for example the translation invariance of the 
coupling ensures that the term arising from the diagonal does not depend on Bt, which 
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is a feature that cannot be reproduced in finite dimension. So while we beheve that 
a theory of double stochastic integrals originating from the variance of a Gaussian 
process could be developed, it is not altogether straightforward and we leave it as a 
future project. 
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